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Abstract 

The southern oscillation index (SOI) is a characteristic of the El Nino phe- 
nomenon. SOI monthly averaged data is analyzed for the time interval 1866- 
2000. The tail of the cumulative distribution of the fluctuations of SOI signal 
is studied in order to characterize the amplitude scaling of the fluctuations and 
the occurrence of extreme events. Large fluctuations are more likely to occur 
than the Gaussian distribution would predict. The time scaling of fluctuations 
is studied by applying the energy spectrum and the Detrended Fluctuation 
Analysis (DFA) statistical method. Self-affine properties are found to be 
pertinent to the SOI signal and therefore suggest power law correlations of 
fluctuations of the signal. Antipersistent type of correlations exist for a time 
interval ranging from about 4 months to about 6 years. This leads to favor 
specific physical models for El Nino description. 



El Nino is one of the most fascinating phenomena in meteorology. Its consequences are 
thought to be wide ranged and sometimes catastrophic. There is no fundamental expla- 
nation, even though much work has been performed on the subject.0 One basic function 
characterizing el Nino is the so-called southern oscillation index (SOI). The cyclic warming 
and cooling of the eastern and central regions in the Pacific Ocean are coupled to distinctive 
sea level pressures. The normalized difference between the pressure measured at Darwin and 
the pressure measured at Tahiti are used to generate the so-called SOI number. When the 
number is negative there is an El- Nino (or ocean warming), but when the number is positive, 
there is a La-Nina effect (or ocean cooling). Indeed the central and eastern Pacific regions 
of the ocean are normally colder than its equatorial location would suggest, mainly due to 
the influence of northeasternly trade winds, a cold ocean current flowing up the coast of 
Chile, and to the up-welling of cold deep water off the coast of Peru. At times, the influence 
of these cold water sources wane, causing the surface of the eastern and central Pacific to 
warm up under the tropical sun. 

There is already a large number of SOI analyses. Here we analyze the southern oscil- 
lation index monthly averaged data for the time interval 1866-2000, i.e. 1612 data points.i'B 
Even though this is a limited number of data points, it is possible to study the amplitude 
and time scaling of the fluctuations of the SOI data by calculating the power spectrum and 
applying the Detrended Fluctuation Analysis statistical method (DFA)! with acceptable 
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error bars. The latter technique has already demonstrated its usefulness in turbulencejl 
in biology for sorting out coding and non-coding sequence behavior in DNA,0 in financial 
analysis for characterizing e.g. currency exchange rates, and in other self-organizing 
critical systems^ for time series with similar lengths. 

We do not aim at forecasting^ here. Forecasting models already exist for El-Nmo events, 
through the sea-surface temperature canonical correlation analysis (CCA) modeljlllll the 
coupled ocean/atmosphere model,0 and the linear inverse model@ to mention a few of them. 

Note also, an excellent, detailed intercomparison of three different forecasting models 
that have predicted moderate warming of the equatorial Pacific sea surface temperatures 
(SSTs) to begin in late summer to early fall of 1986.0 One of the models includes sea level 
pressure data, a quantity that is subject to our study. This is a statistical model, that is 
designed to predict sea surface temperature from prior variations of the sea level pressure. 
Authors draw a conclusion that the key to the success of the model is recognizing the large- 
scale, low-frequency changes in the tropical ocean-atmospheric system. They based their 
conclusion on the estimate of the relative importance of the prior data, included in the 
model, for the subsequent predictive skill of the model. 

The present work stresses the implications of correlations in the amplitude and time 
scaling of the fluctuations of SOI. In so doing, a selected choice between proposed models 
can be envisaged. 



I. DATA 

Data for years 1866-2000 were obtained from Ref. ^| and Ref. |7]. They are plotted as 
function of time (in month units) in Fig. 1. The data consists of 1612 data points and 
represent monthly averaged normalized difference between the pressure measured at Darwin 
and the sea level pressure measured at Tahiti. For the years before 1866 daily measurements 
of the sea level pressure at both stations have been reported to exist and monthly values of 
the southern oscillation index have been calculated!! back to 1841. However, there are gaps 
of a couple of years in the record whence not suitable for our analysis. 



II. DISTRIBUTION OF FLUCTUATIONS 

First, we calculate the fluctuations y i+ i — yi of the SOI signal y(t), group them into bins 
with size equal to the difference between the maximum and the minimum of the fluctuations 
and divide by one hundred, i.e. 0.082. Then we count the number of entries inside each bin. 
The result is the histogram shown in the inset of Fig. 2. The distribution of fluctuations of 
the 5*0/ signal is symmetrical and not Gaussian. 

Next we focus on the tails as they characterize how likely the extreme events are to 
occur. We calculate the empirical probability to observe fluctuations with an amplitude 
larger than some value, Prob(X > x), where x = \yi+i — y%\- The result is plotted in Fig. 
2 (dots). The tail of this distribution is consistent with a power law Prob(X > x) ~ 
showing the reduction of the probability for increasing intensity of the fluctuations. A linear 
least-squares fit yields an estimate /i = 3.30 ± 0.06 when the amplitudes of fluctuations 
are between 1.5 and 2.8. It is difficult to have good estimate for larger amplitudes of the 
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fluctuations (between 3 and 4) due to insufficient statistics. Note, that the //-value is well 
outside the range for stable Levy distributions, < \i < 2, and in particular outside the 
Gaussian distribution (/t = 2). Therefore, large fluctuations are more likely to occur than 
the Gaussian distribution would predict. 

The cumulative probability distribution provides an estimate for the intensity structure, 
e.g. the amplitude scaling of the fluctuations of the 5*0/ signal. The asymptotic power 
law behavior of the intensity distribution of the fluctuations appears to be a particularly 
suitable description of the occurrence of extreme events. In the next section we will focus 
on the time scaling and correlations of the fluctuations of the SOI signal. 



III. SOI TEMPORAL CORRELATIONS 

We apply two methods to study the temporal correlations of fluctuations in the 5*0/ 
signal, spectral and the detrended fluctuation analysis methods. 

The spectral analysis involves calculating of the power spectrum as a Fourier transform 
of the data. To quantify the correlations of the fluctuations in a signal the scaling properties 
of its power spectrum are tested. Assuming a power-law behavior of the energy spectrum 
S(f) ~ f' 13 the self-affine properties of the signal are characterized by the /3-value. 

In Fig. 3 results from the spectral analysis of the SOI data are plotted. For the frequency 
range from 1/5 to 1/64 months -1 the spectrum is consistent with a power-law behavior 
with a spectral exponent (3 = 1.32 ± 0.14. However such Fourier transform analyses fall 
short of precisely showing crossover regimes. To better estimate the crossover and to test 
the correlations using a different approach, we analyze the 5*0/ signal applying the DFA 
technique.! 

The detrended fluctuation function F(r) is calculated following 

1 ( fc+1 ) T N 
F 2 (r) = - £ [y(t) - z(t)} 2 , k = 0, 1, 2, ■■-,(-- 1), (1) 

T t=kr+1 T 

where z(t) = at + b is a linear least-square fit to the data points in a box containing r points. 

The behavior of the averaged though other forms for z(t) can be used0 F 2 (t) over the 
N/ r intervals with length r is expected to be a power law 

<F 2 (r)> ~ r 2 «. (2) 

An exponent « ^ 1/2 in a certain range of r values implies the existence of long-range 
correlations in that time interval as, for example, in fractional Brownian motionEJi^. A 
value of a < 0.5 indicates antipersistence and a > 0.5 implies persistence of correlations. 
The classical random walk (Brownian motion) is such that a = 1/2. 

In Fig. 4, a log- log plot of the function J (F 2 (r)) is shown for the SOI data in Fig. 1. 
This function is close to a power law with an exponent a± = 0.25 ± 0.01 holding for the 
interval time ranging from about 4 to 70 months. In contrast to the Fourier spectrum analysis 
the crossover of the fluctuations function < F 2 (t) > is well defined at 70 months. For time 
scales longer than 70 months, i.e. about 6 years, a crossover to noise-like a<i = 0.05 ± 0.009 
is observed. This suggests antipersistence of the correlations in the fluctuations of the sea 
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level pressure for time lags less than 6 years. Antipersistence of the fluctuations implies that 
a positive fluctuation in the past is more likely to be followed by a negative fluctuation in 
the future. 

Note that a = Hu, Hu being the so called Hurst exponent. The Hurst exponent of a 
signal was first defined in the "rescaled range (R/S) analysis" (of HurstBEl) to estimate the 
correlations in the Nile floodings and droughts. The relationship 

13 = 2Hu + 1 (3) 

has been theoretically proven by Flandrin for fractional Brownian walks.0 The ai\ value for 
the SOI signal is consistent within the error bars with the spectral exponent satisfying the 
above equation. 



IV. CONCLUSIONS 

We have studied the amplitude, time scaling and correlations of the fluctuations of the 
Southern Oscillation Index (SOI). The tail of the cumulative distribution of the SOI 
fluctuations is found to scale with an exponent /i = 3.3 for amplitudes of the fluctuations 
between 1.5 and 2.8, describing the occurrence of extreme events. The energy spectrum 
of the SOI is consistent with power law with exponent (3 = 1.32 for the frequency range 
between 1/5 and about 1/64 months. Since 1 < (3 < 3 the SOI signal is a self-affine fractal£3 
To estimate more precisely the crossover regime and the type of correlations we have applied 
the DFA method. Using the DFA method we find an antipersistent type of correlations for 
time lags less than 70 months. The a-exponent that characterizes the scaling is consistent 
within the error bars with the spectrum scaling. This leads to favor specific physical models 



for El Nino description, as that in Ref. 15 



Our analysis shows that long-range correlations exist between the fluctuations of the 



SOI, e.g. sea level pressure. This supports one conclusion in Ref.[ [L5| pertinent to our 
analysis, i.e. the large-scale low-frequency variations of the global SLP field are responsible 
for the predictive skill of the sea level pressure forced statistical model. The present work 
indicates a hint why it might be so: because of the correlations between the fluctuations. 
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FIG. 1. The monthly values of the Southern Oscillation Index for the period 1866 to April 2000 
representing the standardized pressure difference between Tahiti and Darwin stations (1612 data 
points). 



10 



10' 



A 
O 



10' 



80 
60 

>■ 
o 

CD 

=i 40 

D" 

a> 

LL 

20 




10 
10 




Prob(X>x) ~ x" 



H=3.30 ± 0.06 



10 



x=|y i+1 -y;l 



10 



10 



FIG. 2. The empirical probability to observe fluctuations with amplitude larger than some 
value x = \yi+i — yi\- An asymptotic power law scaling is found for fluctuation amplitudes between 
1.5 and 2.8. Histogram of the fluctuations of the SOI signal is shown in the inset. 
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FIG. 3. The energy spectrum of the SOI data (from Fig. 1). A spectral exponent 
(3 = 1.32 ± 0.14 characterizes the correlations of fluctuations in the frequency range from about 
1/5 to about 1/64 months" 1 . 
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FIG. 4. The DFA-function < F 2 (r) > 1/2 in log-log plot for the SOI data from Fig. 1. Two 
scaling regimes are observed, a\ = 0.25 ± 0.01 and noise-like «2 = 0.05 ± 0.009 with a crossover at 
70 months. 
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